Influence of Supercurrents on Low- Temperature 
Thermopower in Mesoscopic N/S Structures 

J. Zou^, I. Sosnin^, P. Virtanen^, M. Meschke^, 
V. T. Petrashovi, and T. T. Heikkila^ 

^Physics Department, Royal Holloway University of London, 
Egham, Surrey, TW20 OEX, U.K. 
^Low Temperature Laboratory, Helsinki University of Technology, P.O. Box 2200, 

FIN- 020 15 HUT. Finland 



The thermopower of mesoscopic normal metal/superconductor structures has 
been measured at low temperatures. Effect of supercurrent present in nor- 
mal part of the structure was studied in two cases: when it was created by 
applied external magnetic field and when it was applied directly using extra 
superconducting electrodes. Temperature and magnetic field dependencies of 
thermopower are compared to the numerical simulations based on the quasi- 
classical theory of the superconducting proximity effect. 
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1. Introduction 

Thermoelectric properties of hybrid normal- metal/super conductor (N/S) 
structures are strongly modified by the superconducting proximity effect (see 
Befs. [T|2l for review). It was predicted that at low temperatures the ther- 
mopower in mesoscopic N/S structures can be as much as 1000 times larger 
than that in normal metalsP Moreover, the thermopower acquires a phase- 
coherent part as was first observed in Ref-H where thermopower oscillations 
as a function of superconducting phase difference (p were recorded in a geom- 
etry of an Andreev intereferometer (a superconducting loop connected to a 
normal part, see review^^for details). According to theory,!^ this oscillating 
part of thermopower should be antisymmetric in (p and have reentrance in 
temperature dependence with a maximum at the Thouless energy, similar 



J. Zou et al. 



to that of magnetoresistance. The value of thermopower at the maximum 
can reach a few ^V/K. It is worth mentioning that in quantitative mea- 
surements of thermopowerpl where local thermometry technique based on 
proximity resistanc^ was used, the value of only lOOnV/K was reported, 
considerably less than that predicted by theory. While the above qualitative 
behavior has been observed in several experiments, some deviations from it 
have been reported as well. In particular, large thermopower symmetric in 
(/> has been observed,'^^ which is not explained by theory. Close to the su- 
percondicting transition of the superconducting parts of hybrid mesoscopic 
NS junctions, a thermopower due to quasiparticle thermoelectric currents 
has been reported.^^ However, the value of the thermopower observed 
near in these experiments was considerably larger than that predicted by 
theory.^ In structures consisting of a normal ring with only one supercon- 
ducting contact a small thermopower periodic and symmetric with respect 
to magnetic flux through the normal ring has been reported.^ Recently, 
thermal conductivity of Andreev interferometers has been studied as well, 
see Refs. [TKITTn 

In this paper we concentrate on the effect of supercurrents flowing 
through the normal part of the Andreev interferometer, which were first 
highlighted in Ref. [7| We designed our structure so that the supercurrent 
can be created by an applied magnetic flux through the superconducting loop 
or directly from an external power supply using extra superconducting con- 
tacts. The results are compared with the predictions from the quasiclassical 
theory. 



2. Sample fabrication and Measurement 

The samples were fabricated using two-stage e-beam lithography. First, 
a normal part was made of a thermally evaporated 30 nm thick Ag film. 
Then the second layer was made of 55nm thick Al film used as a supercon- 
ductor. To obtain clean interfaces between the layers, the contact area was 
Ar~^ plasma etched before the deposition of the second layer. This in — situ 
etching process produces interface resistance less than 1 for contact area 
of lOOnm X 200nm. Figure 1 shows the geometry of one of the measured 
samples. The sample consists of an ^T-shaped N wire connected to a super- 
conducting loop with two S contacts. Si and 5*2, and to superconducting 
contacts S3 and ^4. On the left side, the N structure is connected to a 

^Although, in the previous case the Thouless energy Et ~ HD/L^ is determined from 
the distance L between the superconductors, and in the latter case from the distance of 
the superconductors to the normal terminals. 
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Fig. 1. SEM micrograph of the measured sample. 



small normal reservoir (which we will call a quasi-reservoir) , which in turn is 
connected to two superconductor electrodes Hi and i?2- By passing current 
from Hi to H2 we were able to vary the temperature of the quasi-reservoir. 
On the right side the N structure is connected to a normal reservoir, which 
is in good thermal contact with massive Au pads so that its temperature is 
fixed by the substrate. 

Measurements were performed in a He'^ cryostat at temperatures from 
0.28 K to 1.5 K with a magnetic field up to 1 mT applied perpendicular to 
the substrate. Resistivity, p, of the Ag film was about 2 /if^cm with diff'usion 
constants, D, about 130 cm^/s. 

Magnetoresistance measurements were performed using conventional ac 
bridge technique. Thermopower measurements have been done using two 
diff'erent methods. In the first method, a heating current, 1^, was a sum of 
dc and 0.5fiA ac currents and it was applied between Hi and H2. Thermo- 
voltage, Vth, between S2 and A''i was measured using a lock- in amplifier on 
the frequency of the ac signal. In the second method, only ac current at the 
frequency / was applied to the heater but the signal was measured by the 
lock- in amplifier at the frequency 2/. A superconducting magnet was used 
to sweep a magnetic field. Zero of magnetic field on all graphs corresponds 
to zero current through the magnet. Since there was a small shift due to 
magnetic field of the Earth, the relative phase of thermopower and magne- 
toresistance oscillations was double checked by repeated measurements to 
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ensure they were measured with the same reference point. 



3. Theory 

Nonequilibrium electrical properties of diffusive superconductor - normal- 
metal heterostructures can be described with the quasiclassical Keldysh- 
Usadel theoryPJ^^I^^^ We model the experimental sample by a system of 
quasi-one-dimensional wires connected to each other at nodes, see Fig. |21 
The primary object to be described is the Keldysh Green's function, 

G(^,x) = (^""(^'"^ , G^ = -UG'')'h, G^ = G^h-hG\ 

(1) 

where G^, G"^, and G^ are the Retarded, Advanced and Keldysh Green 
functions, h = fL^ + fxT-i describes two degrees of freedom of the distri- 
bution function, and fj is the i:the Pauli matrix. Each of these is a 2 x 2 
matrix in the Nambu particle-hole space. In addition, G^ and G"^ satisfy 
the normalization condition (G^)"^ = (G^)^ = I, where / is the identity 
matrix. In this notation, symbols with a "check" (such as G) are chosen to 
represent matrices in the Keldysh space, and symbols with a "hat" (G) are 
matrices in the Nambu space. 

Usadel equation for G{E,x), combined with the boundary conditions 
in the reservoirs and the nodal conditions, is essentially a circuit theory for 
matrix currents 

ii = cjiAiGVG 

flowing in each wire i. Here cTj = e^vp^iDi is the normal-state conductivity, 
Ai the cross section, vp^i the density of states and Di the diffusion constant 
of wire i. The gauge-invariant gradient is denoted as VG = VG — ieA.[%, G], 
where A is the vector potential. Not all parts of this current are conserved, 
but there are leakage terms due to the finite energy and due to inelastic 
scattering, 

V_-ii = e^UF[-iEn + t,,„G]. (2) 

Here fa = fa (8) /. This form is valid in a normal metal where the super- 
conducting order parameter A vanishes. The first leakage term describes 
the dissociation of Cooper pairs entering a normal metal, and it gives the 
finite penetration depth = yhD/E for the proximity effect. The second 
term Sin describes inelastic scattering due to, for example, electron-electron 
or electron-phonon interaction. In what follows, we assume that the latter 
term is finite only in the electrodes and vanishes inside the wires. 
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The matrix currents flowing at different wires are connected via the 
nodal conditions: the sum of the currents /j flowing into each node vanishes 
and the functions G are continuous across each node. 

For the numerics, it is convenient to parametrize the Retarded and 
Advanced Green's function by the complex parameters 9 and such that 

G^ = cosh(0)'r3 + sinh(^)(cos(x)iT2 + sin(x)f'ri) (3a) 
G^ = — cosh(^)f3 — sinh(^)(cos(x)if2 + sin(x)i'ri). (3b) 

Here 9 and x ^^re the complex conjugates of 9 and x- With this parametriza- 
tion, Eq. @ can be broken into four scalar equations, two for the Re- 
tarded/Advanced parts, 

DV'^9 = -2i{E + iO+) sinh(6') + sinh(26l) , (4a) 

DV-jE = 0, DjE = -smh'^{9)vs, vs ^ D{Vx - 2eA/h) . (4b) 

and two for the Keldysh part, 

DV-jL = 0, jL = -VLVfL-7VfT + jsfT, (5a) 

DV-jT = 0, jT = -VTVfT + T^fL + jsfL. (5b) 

Equations (jlj) describe the spectral properties of the system and from their 
solution, one finds the local density of states, energy and position dependent 
diffusion constants, and the spectral supercurrent. The latter two, Eqs. © 
are the kinetic equations, which describe the behavior of the distribution 
functions fi and /t- They can be expressed in forms of static continuity 
equations for the spectral (energy dependent) energy and charge currents jl 
and jT- The coefficients of these equations are obtained from the solutions 
of the spectral equations, 

Vl = ^Tr[l-G^G'^] = ^(1 + [cosh^p - | sinh^p cosh(2Im[x])), (6a) 

Vt = ^Tr[l-G^f3G'^f3] = ^(1 + |cosh6'p + | sinh^p cosh(2Im[x])), (6b) 

T= ^TilG^G^h] = i|sinh0psinh(2Im[x])), (6c) 

js = ^Tr[(G«VG^ - G^VG'^)f3] = Im[- s\n\i\9)vs]/D = Imfo]. (6d) 

Here Vl and Vt are the spectral energy and charge diffusion constants, T 
is an anomalous kinetic coefficient which is finite only in the presence of the 
supercurrent, and js is the spectral supercurrent. 
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Finally, in the reservoirs the Green functions tend into the bulk functions 
of those reservoirs. For Retarded/ Advanced parts in the absence of magnetic 
field this means = —G^ = fa or = in a normal metal and they 
tend to G^ = gf-^ + f {cos{(/))if2 + sin((/))ifi)in a superconductor with the 
superconducting order parameter A = jAle**^, with 

_ 1^1 . _ Asgn£; 

^~ y/{E + iO+)^ - |Aj2' ^y{E + iO+)2 - |A|2 ' 

or 9 = artanh(|A|/E) and x = fp- 

The reservoir values for the distribution functions /l/t given by 

where eV is the potential and T is the temperature of the reservoir. There 
is one exception: for energies below the superconducting gap, Andreev re- 
fiection forbids the energy current into a superconductor. Therefore, there 
the Dirichlet boundary condition of /l = is changed into the vanishing 
of the energy current, ji = into all superconductors. 

If one would like to describe the behavior of the superconducting order 
parameter A in the superconducting parts of the structure, a self-consistency 
equation connecting A with the solution {0,Xi Jl-, fx} could be applied. 
However, in what follows we assume that all the superconductors are reser- 
voirs, such that A obtains its bulk value quickly near the NS boundary. 
Such an assumption works fairly well for our system, but it brings some 
inaccuracy to the exact position of the NS interface as the inverse proximity 
effect suppresses A and 9 close to the interface. Generally these quantities 
have a healing length of the order of the superconducting coherence length, 
^0 = xAW(2A). 

Another consequence of assuming that the superconductors are reser- 
voirs is a boundary condition /^^ = at the normal-superconducting inter- 
face. For energies below the gap, /t decays into the superconductors within 
^0) but above the gap, the decay length is much longer, of the order of the 
charge relaxation length inside the superconductors. 2- The validity of this 
assumption may thus be questionable at high temperatures, where thermal 
quasiparticles inside superconductors play a role. To simulate the effect of a 
very long charge relaxation length, we replaced the Dirichlet boundary con- 
dition for the distribution functions at > A with a Neumann condition 
dxfh = dxfT = 0. This then takes into account the charge-imbalance volt- 
age due to the quasiparticle current entering the superconductor, assuming 
the detailed form of charge relaxation can be neglected. Compared to the 
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{(l^z - 01oop)/2 {(jjz + <?^loop)/2 

Fig. 2. Schematic system considered in the numerics. Major features of the 
experiments can be attributed to the effect of phase gradients across the 
superconducting interfaces on the various observables of the system. The 
phase (pz is fixed by current conservation, Ig = —It; for a left-right symmetric 
system (pz = 0. The phase i;^ioop can be controUed by the flux through 
the superconducting loop, and (ps by the applied supercurrent between the 
contacts S3 and S4. 

Dirichlet condition, this condition results in our geometry to a slightly re- 
duced thermopower at high temperatures, where thermal quasiparticles play 
a role. 

Equations (0]) and © together with the nodal and boundary condi- 
tions need in general to be solved numerically. When the solutions to them 
are found, one can obtain, for example, the observable energy and charge 
currents flowing in wire i from 

lQ,i = [ dEjL,i (7) 

6 JO 

ic,i=^ rdEjT,i. (8) 

6 JO 

The current conservation and Kirchoff laws for these currents follow naturally 
from the corresponding laws for Jlit- 

The numerical results presented in the remainder of this text have been 
obtained by solving the above equations without further approximations. 
We model the experimental setup with the schematic system presented in 
Fig. 121 where seven quasi-one-dimensional wires are connected to each other 
and to two normal-metal and four superconducting reservoirs. The lengths 
and areas of the different wires are estimated from the SEM picture and 
the resistances of each wire, as most of the latter can be separately mea- 
sured. Thus, we estimate Li = 1.44//m, L2 = 1/xm, L3 f» L4 = 340nm, 
and L5 = 460nm, wi = 120nm, W2 = 160nm, 1^3 ~ ti;4 = lOOnm, and 
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= lOOnm, where Lj are the lengths and wi are the widths of the different 
wires. Thickness of the wires was assumed constant. This corresponds to 
the resistances Ri = R2 = A.lVt, ^ = 2.20, and R5 = 30. 

In we included also the resistance 0.8O of the quasireservoir. As the 
parameters for the wires 6 and 7 could not be separately measured, we as- 
sume the system symmetric with respect to the inversion around the around 
the axis lying along wires 1, 5, and 2. In all the numerical curves, we use 
the same set of parameters, and fit only the diffusion constant D and the 
self-inductance L of the superconducting loop, see below. 

We aim to calculate two separately measured observables for the same 
measured structure: the critical supercurrent Ic{T) between the supercon- 
ductors S3 and S4 and the thermopower Q = (V2 — Vs)/iT2 — 7i)|/j=/2=o 
between reservoir 2 and the superconductors. 

As shown for example in Refs. I7|2[ the presence of the spectral super- 
current and the anomalous kinetic coefficient in Eqs. © may lead to a finite 
thermopower. In the limit where the energy scales of the system are far be- 
low A, there is also an approximate relation between the induced potentials 
and the supercurrent flowing in the system, 

^ e i^5(2i?4/3 + R^)Rz/i{Is{Ti) - Is{T2)) 

= 0.58Oe(/5(Ti)-/5(r2)). 

The numerical value in Eq. (13) is specific to the geometry of the present 
system. 

It turns out that the phase difference (p between the two NS interfaces Si 
and S2 and the total flux through the superconducting loop are not linearly 
proportional to each other, but one has to take into account both the kinetic 
and geometric inductance of the loop to find their relation, as was shown in 
Ref . 1251 The former is found using the consistency equation (see for example 
R.ef . I77|) 

h{<t^)=hoop{(t>)- (10) 

To find (/>(<!'), we can as a first step assume that the phase gradient in the 
loop is small and approximate 

7r|A| |A| Aavs Aavs 27rn - (/> + 2^$/^>o 

/loopW' ~ ^^tanh— — — — , ~ ^ , 11 

2 2kBT D D eRN,ioop 

where 27rn — + <I>/<I>o (n S Z) is the gauge-invariant phase difference over 
the superconducting part of the loop, $ the total flux through it and $0 = 
h/(2e) is the fiux quantum. The prefactor in front of vs is found by solving 
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the Usadel equation inside a bulk superconductor. Applying Eq. (|1U|) now 
implies that 



^ 2e , hRiy loon 



tanh 



2kBT 



(12) 



In our structure, assuming Rn,Ioop ~ 20^2, and jA|//cB ~ 2K yields a kinetic 
inductance Lk ^ 20pH at ksT <C |A|. However, one should also take into 
account the geometric self-inductance Lgeom of the SQUID, which modifies 
the relation between $ and the external flux. Hence, the true phase difference 
between the superconductors is found from 

27r^ (mod27r) - </) = "^{Lk + Lgeom)/(0) = '^LI{4>). (13) 

In what follows, for each external flux we calculate the phase (j) by 
solving Eq. ((T^ numerically. It turns out (see below) that the results are 
well fitted with the loop inductance L « 50 pH. For a circular loop of radius 
R ~ 2.4 /im and cross section Ai = 35nm x 400 nm, we would estimate a 
geometric self-inductance -Lgeom = HoR\^'n.{SR/ Ai) — 7/4] = 21 pH, so that 
the sum Lk + Lgeom = 40 pH is comparable to the fitted value. 

It is straightforward to verify that the Usadel equations (|3HS1) remain 
invariant under the transformation 

-X, A -A, /t ^ -/t, j't ^ -Jt , (14) 

also inside the superconductors. Moreover, it turns out that this symme- 
try is shared by Eq. (fT^ and the self-consistency condition, provided that 
A = jAje*"^ 1-^ lAle"*"^. This implies that for any given solution of the prob- 
lem, there is a second solution with inverted magnetic and electric fields, 
corresponding to charge and supercurrents flowing in the opposite direction. 
For the thermopower this implies that between any pair of reservoirs. 



AT(0) 



_ -AV{- 



AV{- 



-ia=o 



(15) 

lc=0 



i.e., that its phase oscillations are always antisymmetric. This is a general 
symmetry of the quasiclassical model in a static situation, valid in a flnite 
magnetic field and also when superconductors are treated self-consistently. 

This symmetry is in contrast with the cos((/))-dependent thermopower 
measured in Refs. I4pi()[ or the constant offset thermopower measured in 
Ref. 1141 At present we do not know a reason for this discrepancy. 
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Fig. 3. Critical current Ic of the SNS junction as a function of temperature. 
Dots: measured Ic-, solid line: fit to the theory (see text). Inset: Differential 
I — V curve of the SNS junction at T = 0.28 K measured using S^-S^. 

4. Experimental data 

Below we present experimental data for thermopower and resistance 
measurements in the presence of temperature gradients created by heater 
currents. To avoid heat leak into electrodes Hi and H2 we kept Ih < Ic, 
where Ic is the critical current of the superconducting transition for heater 
electrodes Hi and H2- We measured current- voltage characteristic (I-V) of 
H1-H2 and found Ic = 25^A, so Ih smaller than lOfiA was used to ensure 
contacts Hi and H2 are in the superconducting state. 

To characterize the SNS junction, we measured its critical current at 
different temperatures. Inset of Fig. 01 shows differential I — V curve of the 
SNS junction measured using contacts 53-54 at r=0.28K. The jump at high 
currents corresponds to the superconducting transition of the parts of Al 
contacts 53-^4, while the second jump at smaller currents corresponds to 
the superconducting transition of normal parts L3, L5, and L4. The critical 
current was measured as the current value at dV/dI=Rj\f /2, where Rn ^ 
is the value of resistance of L^-L^-L^i in normal state. Asymmetry in 
dV/dl curve can be attributed to the difference in Joule heat released in 
the wires when the superconducting transition is approached from normal 
state (negative currents in Fig. 3, inset) and that when it is approached 
from the superconducting state (positive currents in Fig. 3, inset). The 
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temperature dependence of Ic is plotted in Fig. El A numerical best fit, 
based on the quasiclassical theory presented above, to this data is given 
as a solid line, which also includes the temperature dependence of the gap 
A(r). We used low temperature value of the gap A(0)//cb ~ 2K estimated 
from the superconducting transition temperature of Al wires. In this fit, 
we found the Thouless energy Et ~ SOmK/A;^ that best corresponds to the 
measured temperature dependence. The fitted Thouless energy corresponds 
to a diffusion constant D « 85cm^/s, in accordance to the value estimated 
from resistivity. 




Fig. 4. (Color online): Amplitude of magnetoresistance oscillations measured 
with current probes Hi-Ni and potential probes H2-N2. Circles correspond 
to resistance minimum, squares, to maximum. 



Figure shows temperature dependence of the resistance measured us- 
ing current contacts Hi-Ni and potential contacts H2-N2 (see Fig. ^ at 
different values of magnetic flux through the S loop. The bottom curve is 
for $ = and the top one is for ^ = ^o/2. At a fixed temperature the 
resistance oscillates between these two values as a function of the applied 
magnetic flux, see Fig. 5b. The amplitude of the oscillations is affected by 
the screening effect discussed in Ref. 25 261 The screening effect was also seen 
directly as deviation of the shape of magnetoresistance oscillations from sinu- 
soidal form, see Fig. |SI Figure IHl shows magnetoresistance and thermopower 
oscillations at the same temperature. The magnetoresistance was measured 
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using current probes Hi-Ni and voltage probes H2-N2. Thermovoltage was 
measured using heater current applied to heaters H\-H2. For thermopower 
measurements presented here the heater current was a sum of 4 /xA dc and 
0.5 /xA ac components. The signal was then measured by the lock-in ampli- 
fier at the frequency of ac modulation between contacts S2-N1. Note that 
thermovoltage oscillations are antisymmetric with respect to the direction 
of magnetic field, in contrast to magnetoresistance oscillations (this is seen 
directly from phase shift of 7r/2 for thermovoltage oscillations compared to 
that of magnetoresistance). 




H(G) 

Fig. 5. Top: Thermovoltage oscillations measured at T = 0.28 K with 
heater current H1-H2 and voltage probes S2-N1. Bottom: magnetoresistance 
oscillations measured at T = 0.28 K using current probes Hi-Ni and voltage 
probes H2-N2. 

Non-sinusoidal shape of magnetoresistance oscillations indicates the pres- 
ence of supercurrent in the normal part of the hybrid superconducting loop. 
Indeed, direct measurement of I-V curves of the central A^-wire using con- 
tacts S^-Si confirms this (see Fig. ^ and gives us an estimate of the mag- 
nitude of this supercurrent at different temperatures: the SNS structures 
formed by the superconducting loop and its contacts to the normal-metal 
wires are very similar to that formed by the S^-S^ contacts. We may thus 
estimate the effect of screening from the shape of the thermovoltage oscil- 
lations, using the fitted value of the supercurrent. Figure IHl shows the mea- 
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Fig. 6. (Color online): Thermovoltage oscillations (filled circles) measured 
at T = 0.28 K with heater current H1-H2 and voltage probes S2-N1. The- 
oretical fits to the shape of the oscillations with different values of the loop 
self-inductance: L = (solid line), L = 50 pH (dashed line) and L = 100 
pH (dash-dotted line). 

sured thermovoltage oscillations along with three theory curves calculated 
with different loop inductances. Note that the amplitude of the theory curves 
are scaled to the experimental data, so one should only pay attention to the 
shape of the oscillations. The temperature dependence of the measured ther- 
mopower amplitude is compared to the theory in Fig. EH — the inductance 
fit is mostly sensitive to the S loop screening parameter f3 = LIc/^o, inde- 
pendent of the thermopower amplitude. We find the best fit with L ~ 50pH. 
This is close to the lower limit of L ~ 40pH estimated from the shift of the 
magnetoresistance oscillations in the presence of the supercurrent applied 
from a power supply connected to 83-84.'^ There was no hysteresis in either 
magnetoresistance or thermopower oscillations versus applied magnetic flux. 
For L=50pH and /c=6^A the screening parameter /? = 0.15. According to 
theory,!^ one sees hysteresis when /3 > l/vr, so the above estimation of L is 
consistent with the absence of hysteresis. 

In Fig. we plot the measured thermovoltage oscillations at three dif- 
ferent temperatures, along with similar normalized theory curves as in Fig. El 

■^Note that this only gives a lower limit, as only part of the supercurrent between S3 
and Sa enters the superconducting loop. 
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Fig. 7. (Color online): Measured thermovoltage oscillations at three different 
temperatures, from bottom to top: T = 0.28 K, T = 0.5 K and T = 0.7 
K. The curves have been offset for clarity and normalized to the values 
at $ = $o/4. The solid lines on top of the data are the corresponding 
theory curves calculating with the inductance L = 50 pH. One should pay 
attention only to the shape of the curves. The temperature dependence of 
the measured thermopower amplitude is separately compared to the theory 
in Fig. Uni 



with L = 50pH. The temperature dependence of the oscillation amplitude is 
compared to the theory separately, see below. 

Figure IHl shows a temperature increase of the central part of interferome- 
ter Tm as a function of dc heater current at different temperatures. This was 
measured by comparing temperature dependence of the resistance of the N 
wire, measured using current contacts S3-N1 and potential contacts S2-N2, 
to the dependence of that on the dc heating current (see also Refs. I8|9|ll|l . 
Measurements using the critical current of an SNS junction give similar de- 
pendences. The solid lines represent an approximation of low heating current 
part of the curve to a dependence Tm = Tq + al^, where a is a constant, 
which is valid when al^ <C Tq. The values of 2a=(PTm/dl'^ are shown in 
Fig. ini However, in order to calculate thermopower we need to use values of 
2a*=(PTh/dl'^ , where is the temperature of the "hot" end of the struc- 
ture, which in our case is the normal quasi-reservoir (see Fig. 1). Numerical 
simulation based on the actual geometry and values of resistances Ri - R5, 
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Fig. 8. (Color online): Temperature rise in the middle of the interferometer 
as function of heater current H1-H2, measured at T = 0.28 K, 0.5 K, 0.7 K, 
and 0.9 K 

which also took into account the error in measuring using the above 
method due to the difference in the interferometer resistance heated uni- 
formly compared to that in the temperature gradient, showed that a* ~ 4a, 
which was accurate within 10% in the whole temperature range used. Now 
we can convert the measured thermovoltage into thermopower. Thermo- 
voltage measured at the frequency / of the ac modulation part of heating 
current is given by 

= -TF^-c = Q-JT^-c = Q'^^*IdJac. (16) 
al di 

where Idc and lac are dc and ac components of the heating current, respec- 
tively. Thermovoltage measured at the frequency 2/ in case of Idc = can 
be presented as (see also Ref. |H1). 
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Note that the extra factor 1/2 in (17) arises from the relation I^^cos^ft = 
(l/2)/^^+(l/2)/^(,cos2/t, so that lock-in amplifier output at 2/ is propor- 
tional to (l/2)/2^. 

We have measured thermopower using both methods and found the two 
in a reasonable agreement. For example, for T = 0.28 K we found Qf = 1-2 
//V/K and Q2/ = 0.9 //V/K. In the range of heating currents Ih < 10/iA, both 
Vf and V2/ always had symmetry of sine/) and showed no phase-independent 
component observed in other experiments. 

Figure shows the value of thermopower as a function of bath tem- 
perature. The thermovoltage in this case was measured at 2/ with I^c = 
and lac = 6//^. Magnetic field during this measurement was such that 
$ = ^0/2. Q was then calculated using (17) and values of 2o found from 
Fig. 1^ by a linear interpolation between the experimental data points. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 

T(K) 

Fig. 9. Temperature dependence of the second derivative of the effective 
temperature rise with respect to heating current. Points are experimental 
data. Lines represent linear interpolation used in calculation of Q. 

The temperature-dependent Q calculated from the theory is presented 
as lines in Fig. IIUI Similar to the experimental data, the theory predicts 
a reentrant thermopower, with the maximum slightly above the Thouless 
energy. This is below the maximum point ~ 0.47 K found in the experi- 
ments. Increasing the Thouless energy in the simulations would lead to an 
improved fit of the peak position (see the dotted curve in Fig. [rU)) . but then 
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the critical current data cannot be understood. In the presence of screening 
(finite L), at the lowest temperatures the phase does not reach the point 
where the maximum thermopower would be obtained, and thus the result- 
ing thermopower is reduced, and the position of the maximum thermopower 
is slightly shifted to the right. However, to get the maximum near the ex- 
perimentally observed value, a much larger value of L would be needed, and 
in this case the calculated Q{T) is much wider than in the experiments. 




Q[ , , , , 

0.2 0.4 0.6 0.8 



T (K) 

Fig. 10. (Color online): Thermopower as a function of temperature. The 
circles are experimental data. Lines are numerically calculated amplitudes, 
for Ex/kB = 0.05 K and different values of the inductance: L = (solid 
line), L = 50 pH (dashed line), and L = 500 pH (dash-dotted line). For 
comparison, similar data calculated for Ex/kB = 0.4 K and L = is plotted 
as a dotted line. 

In the geometry of our experiment the effect of supercurrent can be 
measured directly. Figure ^2 shows the dependence of thermovoltage on 
magnetic field measured at different values of dc current between contacts 
53-5*4. At small currents an extra shift in thermovoltage occurs. The precise 
shift depends on the ratio of the supercurrents entering the loop and flowing 
in the central wire. According to the simulations, the previous dominates at 
the temperatures where the measurements were made, and as a result the 
phase shift is close to aicsm{Is / Ic)- 

Measuring Vj as a function of at = showed linear dependence 
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(see Fig. ^jand Eq. (|18jl ). At higher values of current through 53-54 the 
total current in the central N-wire (applied dc current plus screening current 
due to magnetic flux) exceeds the value of critical current of this part 
see Fig. 4). In this dc voltage will appear between thermovoltage 

probes Sl-Nl, which is also temperature dependent. Since temperature is 
modulated at the frequency /, so is the above constant voltage, leading to 
an extra contribution to signal measured by lock-in amplifier. 

Let us concentrate on the linear response to Is- The effect of Is can be 
estimated using Eq. (13). Assuming the supercurrent - phase relation in the 
N wire is not strongly temperature dependent, we have Is{T) = Ic{T)^{(l}), 
where ~ sin(0), Eq. (13) can be rewritten as 

Q{Is) = 0.58f]^^ = 0.58f]^^ sin0 = 0.58J7^^-^. (18) 
^ dT dT ^ dT Ic{T) ^ ' 

Substituting dl^/dT = 8.7 /^A/K at T = 0.28 K found from Fig. Elinto (18) 
we get Q = 2.5 /iV/K for = 3 /iA. Comparing with experimental data we 
have Vf{2>^A) = 20 nV at I^c = 5 /iA and lac = 0.5 /zA, which corresponds to 
Q = 0.7 /iV/K. The discrepancy should be attributed to a more complicated 
dependence between Ig and (j) due to the presence of the 5 loop. 
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Fig. 12. (Color online): Thermovoltage as function of dc supercurrent ap- 
plied between S3-S4 at T=0.28K. Blue: magnetic field corresponds to a zero 
flux through the S loop. Red: magnetic field corresponds to a half flux 
quantum. 

5. Analysis and Discussion 

We have explored the dependence of the thermopower of an Andreev in- 
terferometer on supercurrents created by magnetic field and applied directly 
from external power supply. The experimental thermopower was modelled 
theoretically using the numerical calculations with the actual geometry and 
fitting only the Thouless energy corresponding to the distance between the 
superconductors. We find that the thermovoltage oscillates as a function of 
the magnetic flux through the superconducting loop. The shape of these 
oscillations can be well described by the theory once the loop inductance is 
taken into account. The theory predicts correct order of magnitude for the 
thermopower. However, the calculated temperature dependence does not fit 
experimental data. The discrepancy can be attributed to the experimental 
error in calculating Q due to the ambiguity in determination of true temper- 
ature difference across the interferometer using only one thermometer in the 
middle of the structure as opposed to two separate thermometers at hot and 
cold ends of the interferometer. Another possible reason is a complicated 
geometry with its extra superconducting electrodes 5*3 and 5*4 in addition to 
the S loop, and especially the ambiguity in determining the parameters for 
the sample geometry. 
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The effect of supercurrent has been measured directly using extra super- 
conducting electrodes. The experimental dependence is linear in agreement 
with theory. The absolute value of thermopower due to this supercurrent 
found experimentally is smaller than that predicted by theory. This dicrep- 
ancy may at least partially be attributed to a more complicated relation 
between Q and Ig due to the presence of the S loop. When both externally 
applied supercurrent and screening current due to magnetic flux through the 
S loop are present, the situation becomes more complicated (see red curve 
in Fig. I12|) . In particular, when the total current approaches the critical 
value of the N wire, the dependence of the thermovoltage on Ig becomes 
strongly nonlinear. This cannot be accounted for completely by only the 
above mentioned extra contribution to measured signal. This case will be 
discussed in detail elsewhere. 

When Ih < Ic, the quasi-reservoir is heated locally so that its distri- 
bution function is close to equilibrium. The main heat transport channel is 
through the A'^-wire into reservoir, so that in this regime the temperature 
gradient is well defined. We concentrated on this regime to compare ob- 
tained results with the theoretical calculations. When Ih > Ic (or when the 
heater electrodes made of a normal metal) heater contacts turn normal and 
the quasi-reservoir distribution function has a nonequilibrium form due to a 
long energy relaxation length in mesoscopic conductors. The measurements 
in this regime will be reported elsewhere. 
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